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ABSTRACT 



The South China Sea (SCS), the largest marginal sea in the West Pacific Ocean, is 
separated from adjacent oceans by a chain of islands. The deepest water is confined to a 
bowel-type trench, and the maximum depth is approximately 5,000 m. Most of the existing 
studies on the seasonal and interannual variability have been based only on surface 
temperature data. However a primary need is an understanding of the three-dimensional 
thermohaline circulation. The minimum curvature with spline method was used to establish 
a three-dimensional monthly-varying gridded data from the Navy’s Master Oceanographic 
Observation Data Set (approximate 189,000 profiles), covering the area of 5°N - 25°N and 
105°E - 125°E and iBrom the surface to 400 m depth. For temperature, profiles were binned 
into 204 monthly data sets fi'om 1968 to 1984 (17-ye^s). Because of the paucity of salinity 
data, salinity profiles were binned into 12 climatological monthly data sets, and the monthly 
climatological mean was computed. After the gridded data set had been established, both 
composite analysis and the Empirical Orthogonal Function (EOF) analysis were used to 
identify the major thermohaline features. The first EOF mode accounts for 26.7% of the 
variance and represents the seasonal variation. The second EOF mode accounts for 17.7% 
of the variance and represents the interannual SCS warming/cooling phases. Furthermore, 
the P-vector method was used to invert three-dimensional velocity fields from the analyzed 
temperature and salinity data. Important dynamical processes, including the Kuroshio 
intrusion, the western boundary current (counter-current), the cross basin current (under 
counter-current), the mesoscale eddies, and the basin gyre are identified. 
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I. INTRODUCTION 



The South China Sea (SCS), the largest marginal sea in the West Pacific Ocean, is 
separated fi'om adjacent oceans by a chain of islands. It contains the broad shallow Sunda 
Shelf on the south, a shelf extending from the Gulf of Tonkin to the Taiwan Strait on the 
north, a deep sea basin in the center, and the Dangerous Ground with numerous reef islands 
over the southeast (Figure 1.1). It has a bottom topography that makes it a unique semi- 
enclosed ocean basin that is overlaid by a pronounced monsoon surface wind. Extensive 
continental shelves (less than 100 m deep) are found on the western and southern parts, while 
steep slopes with almost no shelves are found in the eastern part of SCS. The deepest water 
is confined to a bowel-type trench and the maximum depth is approximately 5,000 m. 

The SCS's physical oceanographic condition has been investigated for several decades. 
Studies show that eddies in the SCS are predominantly cyclonic in winter and anti-cyclonic 
in summer, with sizes ranging fi’om small to meso-scale, and are both cold and warm [Huang 
et al., 1994]. A survey by Wyrtki [1961] revealed complex temporal and spatial features of 
the surface currents in both the SCS and the surrounding waters. By April, the northeast 
monsoon has decayed and the main SCS current consists of two large cyclonic eddies in the 
western SCS, one in the north and the other in the south (Figure 1.2). With the beginning of 
the southwest monsoon in May and June, a wind drift current is formed off the coast of 
Vietnam. The westward intensification of this current is clearly visible. The wide uniform 
drift in the northern SCS shows a deflection of the current to the right of the wind. A large 
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part of the water passes south of Taiwan and into the root of the Kuroshio. The southern SCS 
is occupied by an anticyclonic eddy. 

Based on limited data sets, studies show that eddies in the SCS are predominantly 
cyclonic in winter and spring and anti-cyclonic in summer, with sizes ranging from small to 
meso-scale [Huang et al., 1994], Both cold and warm eddies exist in the SCS. Dale [1956] 
and Uda and Nakao [1972] reported a cold eddy off the central Vietnamese Coast in summer, 
and Nitan [1970] found a cold eddy located to the northwest of Luzon. Reports from the 
South China Sea Institute of Oceanology [SCSIO, 1985] indicate that in the central SCS a 
warm-core eddy appears in summer and winter, but it is more closer to Vietnam in summer 
at the surface. In addition, a cold-core eddy was detected in the central SCS from the 29th 
December 1993 to 5th January 1994 according to analysis of the TOPEX/ POSEIDON data 
[Soong, et al., 1995]. 

Recently, Chu et al. [1997a] and Chu and Chang [1997] identified a central SCS 
surface warm-core eddy in mid-May from a more complete historical data set - the U.S. 
Navy's Master Observational Oceanographic Data Set (MOODS). From the composite 
analysis of the U.S. National Centers for Environmental Prediction (NCEP) monthly SST 
fields (1982-1994), Chu et al. [1997b] found that during the spring-to-summer monsoon 
transition (March to May) a warm anomaly (greater than 1.8‘’C) is formed in the central SCS 

t 

at 112°-119°30E, 15°-19° 30N. This warm eddy appears in the central SCS (west of Luzon 
Island) during the late spring season and strengthens until the onset of the summer monsoon 
(mid-May) and then weakens and disappears at the end of May. Although its size and 
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intensity varies, the SCS warm pool releases large moisture and heat fluxes into the 
atmosphere and in turn affects the monsoon circulation. 

Most of the existing studies on the seasonal eddy variability were based on the surface 
data. What is the three-dimensional thermal structure and circulation pattern of the SCS 
eddies? In this thesis, we will use the minimum curvature spline method [Wessel and Smith, 
1992] to establish a three-dimensional monthly varying gridded data from the MOODS, 
covering the area of 5°N - 25°N and 105°E - 125°E and from the surface to 400 m depth. 
For temperature, we binned the profiles into 204 monthly data sets from 1968 to 1984 (17 
years). For salinity, a monthly climatological analysis was apphed because of the paucity of 
data. After the gridded data set had been established, we utilized both composite analysis and 
the Empirical Orthogonal Function (EOF) analysis identifing the major thermal features of the 
mesoscale eddies. Furthermore, three-dimensional velodty fields were inverted from analysis 
of the temperature and sahnity data. 

The outline of this thesis is as follows; Chapter II depicts the South China Sea 
oceanography including the water masses, monsoon effect and mesoscale eddies. Chapter III 
describes the U.S. Navy’s Master Oceanographic Observation Data Set (MOODS). Chapter 
IV describes the establishment of gridded data, and Chapter V depicts the composite analysis 
to obtain tempjerature and salinity seasonal variations. Chapter VI describes the EOF method 
to obtain nonseasonal variabilities. Chapter VII presents the SCS circulation using the P 
vector inverse method, major current systems and seasonal variation, and Chapter VIII 
includes our conclusions. 
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Figure 1.1 Bathymetry of the South China Sea 
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Figure 1.2 Observational Surface Currents of the Southeast Asia Waters in April 

(From Wyrtki, 1961) 
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n. THE SOUTH CHINA SEA OCEANOGRAPHY 



A. WATER MASSES 

On the western side of the Luzon Strait, the SCS water masses can be divided into 
those four types [Wyrtki, 1961]; (1) the warm, less salinity Surface Layer Water, (2) the high 
salinity (34.8-35.2 psu) Subtropical Lower Water, (3) the low salinity (34.1-34.5 psu) 
Intermediate Water, and (4) the original Pacific Deep Water ( 34.65-34.7 psu). The salinity 
maximum of the Subtropical Lower Water is approximately 150 meters deep and the salinity 
minimum of the Intermediate Water is approximately 400 to 500 meters deep. At the Luzon 
Strait these two water masses move in opposite directions and the depth of the inversion is 
about 200 meter. The center of this inversion lies relatively uniformly at a temperature 16°C. 

Wyrtki [1961] also found that during the winter monsoon season surface water 
masses fi'om the Pacific Ocean are transported into the SCS through the Luzon Strait, chiefly 
along the western side of the SCS in depths between 400 and 900 meters. During the summer 
monsoon season these conditions are reversed. 

Waters fi’om the Western Pacific enter the northern SCS through the Luzon Strait and 
after mixing and modification form distinct water masses. Continental fresh water runoff is 
also very significant. Su and Weng [1994] classified eight water masses in the northern SCS. 
They grouped these water masses into three classes: (i) runoff diluted type of Near shore 
Diluted Water Mass (F)\ (ii) the shallow sea modified type composed of Coastal Mixed Water 
Mass (M), Warm Surface Water Mass (WSJ, Surface Water Mass (S) and Surface-Subsurface 
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Mixed Water Mass (5C/); and (iii) the deep sea type, composed of Subsurface Water Mass 
(U), Subsurface-Intermediate Mixed Water Mass (UI) and Intermediate Water Mass (I). 
Distribution of these water masses at the surface are nearly parallel to the coastline (Figure 
2 . 1 ). 

Both water masses WS and S originate from West Pacific Surface Water Mass. 
Various meteorological conditions modify this water mass and ultimately generate the high 
temperature and low salinity WS water mass and the low temperature, low salinity S water 
mass. Mixing modifies the West Pacific Subtropical Subsurface Water Mass to form the SU 
water mass, which is similar to the S water mass. The U water mass, distributed at 
approximately 150 m depth, retains the high salinity layer of the West Pacific water. The Ul 
water mass is the interface between the deep high salinity layer and low salinity surface layer. 
The / water mass is an oceanic type water mass with the least modification from its source 
water mass and characterizes low salinity, distributed down to 700 m. The F water mass is 
a mixture of runoff and Near shore water mass. The M water mass undergoes the largest 
modification, mbdng with SU, WS and S water masses, respectively, in varying seasons. Table 
2. 1 hsts the characteristics of the water masses. 
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Water Mass 


Feb. 


May 


Aug. 


Nov. 


F T(°C) 

Salinity 




25.5-27.8 

<30.0 


29.0-30.5 

<30.0 




M TCC) 

Salinity 


13.2-16.0 

31.20-33.60 


22.5-27.6 

31.16-33.80 


28.9-30.3 

33.03-32.91 


21.3-24.4 

31.32-34.02 


WS TCQ 
Salinity 


22.3-24.4 

34.08-34.67 


26.2-30.0 

33.90-34.69 


28.0-30.1 

33.03-34.40 


25.5-26.1 

33.68-33.93 


S TCQ 

Salinity 


20.4-22.3 

34.43-34.80 


22.2-26.2 

33.85-34.62 


23.5-28.0 

33.65-34.62 


22.5-25.5 

33.97-34.60 


SU TCQ 
Salinity 


17.5-20.4 

34.18-34.75 


19.8-22.2 

34.16-34.71 


21.0-23.5 

34.32-34.70 


19.0-22.5 

34.42-34.76 


U TCQ 

Salinity 


15.5-17.5 

34.28-34.78 


14.5-19.8 

34.44-34.88 


15.2-21.0 

34.35-34.81 


14.5-19.0 

34.50-34.78 


UI TCQ 

Salinity 


11.0-15.5 

34.46-34.60 


10.0-14.5 

34.40-34.70 


11.0-15.2 

34.41-34.79 


11.0-14.5 

34.42-34.62 


I TCQ 
Salinity 




5.3-10.0 

34.33-34.64 


5.6-11.0 

34.34-34.52 


6.4-11.0 

34.30-34.50 



Table 2.1. Ranges of Temperature and Salinity of Water Masses 
in the northern South China Sea (from Su and Weng, 1994) 



B. MONSOON EFFECT ON THE CURRENT SYSTEM 

The entire region of the South China Sea is under the influence of the monsoon 
system, and in the absence of major oceanic inflow, the currents undergo a seasonal reversal 
of direction. 

1. Summer monsoon 

During the summer monsoon season, a warm and weaker southwest wind blows over 
the SCS. The southwest-northeast oriented jet stream is positioned at the central SCS 
(Figure 2.2). Such a typical summer monsoon pattern lasts nearly four months (mid-May to 
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mid-September). The summer monsoon forces a wide, uniform northeasterly flowing surface 
current over most of the basin (Figure 2.3). Westward intensification generates a strong 
current along the western side of the basin, off the coasts of Thailand and Vietnam. The 
larger part of the surface flow passes south of Taiwan and merges with the northerly flowing 
Kuroshio; a smaller amount flows north through the Taiwan Strait. Wyrtki (1961) reported 
the existence of a counter current on the eastern side of the basin, with a southwesterly flow 
along the Borneo coast. This flow turns back northerly in the region of Natuna Island, in the 
southern SCS, forming an anticyclonic gyre. Wyrtki hypothesizes that an insufficient supply 
of water into the basin fi’om the Java Sea produces this gyre. In September, these movements 
decrease and in October the northeast monsoon begins blowing with considerable strength, 
resulting in the complete reversal of the current fields. 

2. Winter Monsobn 

During the winter monsoon season, a cold northeast wind blows over the SCS; as a 
result of the Siberican high pressure system located over the east Asian continent. Radiative 
cooling and persistent cold air advection maintain cold air over the SCS. Tne northeast- 
southwest oriented jet stream is positioned at the central SCS (Figure 2.4). Such a typical 
winter monsoon pattern lasts nearly six months (November to April). The mean northeast 
wind stress over the Western Pacific deflects the water masses of the North Equatorial 
Current to the south, causing them to enter the SCS through the Luzon Strait. Water from 
the East China Sea moves southerly through the Taiwan Strait and enters the SCS. This 
southerly moving water turns southwesterly and flows along the coast of China and Vietnam, 
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with a remarkable westward intensification (Wyrtki, 1961). In the center of the SCS, a 
northerly flowing counter current develops over about 10 degrees of latitude (Figure 2.5). 

During the early period of the northeast monsoon, not all of the southerly flowing 
water can flow into the Java Sea, since the southeast monsoon is still blowing weakly over 
this region. This situation produces a northward flowing counter current along the coast of 
Borneo. The water flows as far north as Luzon, where it turns southward and merges into 
the main current, forming a closed cyclonic gyre. 

During December the northwest monsoon is fully developed and the currents are their 
strongest, exceeding 100 cm/s oflF the coast of Vietnam. The northward flowing counter 
current in the central SCS is located farther to the east, due to the strength of the main 
current. A weak southward flowing coastal counter current forms along the northwest coast 
of Borneo. 

C. MESOSCALE EDDIES 

Seasonal occurrence of the SCS eddies have been reported by several authors. 
Huang and Wang [1994] mentioned that eddies in the SCS are predominantly cyclonic in 
winter and anticyclonic in summer, with sizes fi’om small (100 km) to medium scale (500 
km). Dale [1956] and Uda and Nakao [1972] reported a cold eddy off the central Vietnam 
coast in summer. Nitani [1970] found a cold eddy located northwest of Luzon in summer. 
Reports fi'om the South China Sea Institute of Oceanology (SCSIO) [1985] indicated that a 
warm-core eddy appeared in summer and winter in the central SCS, but in summer, it was 
closer to "Vietnam at the surface. Recently, a cold-core eddy was detected in the central SCS 
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between December 29, 1993, and January 5, 1994, from the analysis of TOPEX/ POSEIDON 
data [Soong et al. 1995], Chu et al. [1997a] and Chu and Chang [1997] identified the 
existence of a central SCS surface warm-core eddy in mid-May from a historical data set: the 
U.S. Navy's Master Observational Oceanographic Data Set (MOODS). From the composite 
analysis of the U.S. National Centers for Environmental Prediction (NCEP), monthly sea 
surface temperature (SST) fields (1982-1994), Chu et al. [1997b] found that during the 
spring-to-summer monsoon transition (March-May) a warm anomaly (greater than 1.8°C) 
was formed in the central SCS at 112°-1 19°30E, 15°-19°30'N. 

D. AN AIR-SEA FEEDBACK SCENARIO 

The circulation at intermediate to upper layers of the SCS are mainly influenced by 
monsoons. The prevailing summer (winter) monsoon results in an cyclonic (anti-cyclonic) 
circulation of the surface layer. The reason may be that the SCS is large, semi-enclosed and 
only slightly affected by the Kuroshio [Huang et al.,1994]. Chu and Chang [1997] 
investigated the historical SST data from the MOODS, and proposed an air-sea feedback 
scenario for the central SCS warm-core eddy formation (Figure 2.6). The anticyclonic surface 
winds prevail in the central SCS from late winter to spring. This anticyclonic wind forcing 
generates downwelling in the central SCS and in turn prevents the cold deep water from being 
adverted to the surface. This will promote the formation of the warm-core eddy in the central 
SCS. After the summer monsoon onset, an atmospheric surface cyclone occupies the central 
SCS, which entrains the deep cold water into the surface mixed layer. This upwelling effect 
finally destroys the central SCS warm-core eddy. 
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Figure 2.1 Distribution of water masses at the sea surface (a-d) and bottom layer 
(e-h) in the northern South China Sea. a, February; b, May; c, August; 
d, November; e, February; f, May; g, August; h, November; 
see text for acronyms (from Su and Weng, 1994) 
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Figure 2.2 Mean Atmospheric Surface Streamline Analysis 
for summer monsoon (From Cheang, 1980) 
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Figure 2.3 Observational Summer (June) surface Circulation of 
the Southeast Asia Waters (From Wyrtki, 1961) 
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Figure 2.4 Mean Atmospheric Surface Streamline Analysis 
for winter monsoon (From Cheang, 1980) 
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Figure 2.5 Observational Winter (December) Surface Circulation of 
the southeast Asian Waters (From Wyrtki, 1961) 



17 



Central SCS Surface Atmospheric Anti-cyclone 



Central SCS Downwelling 



Isolating Deep Cold Water 



Promotes Surface Warm-Core Eddy 



Atmospheric Low Pressure at the Central SCS 



Atmospheric Monsoon Onset over SCS 



Downwelling Effect Disappears 



Figure 2.6 An Air- Sea Feedback Scenario for the SCS Warm-Core Eddy Life Cycle 

(From Chu and Chang 1997) 
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m. MASTER OCEANOGRAPHIC OBSERVATION DATA SET (MOODS) 



The MOODS is a compilation of observed ocean data worldwide consisting of (a) 
temperature-only profiles; (b) both temperature and salinity profiles; (c) sound-speed profiles; 
and (d) surface temperature (drifting buoy). These measurements are, in general, irregular 
in time and space. The data are from the early nineteen century to the present and include 
189,059 temperature profiles (Figure 3.1) and 14,004 salinity profiles (Figure 3.2). In this 
study, we analyze the temperature profiles obtained during 1968-84 from a variety of 
instruments. Our study domain includes the area 5°N to 25°N and 105°E to 125°E. The 
primary editing procedures included removal of profiles with obviously erroneous location, 
profiles with laige spikes, and profiles displaying features that do not match the characteristics 
of surrounding profiles. In shallow water, this procedure can be partially automated, but it 
also involves subjective interpretation because of the under sampling of MOODS, compared 
to the spatial and temporal variability of the ocean. 

The temporal and spatial distribution of the MOODS data are irregular. The temporal 
distribution (Figure 3.3) indicates that most observations were obtained during the Vietnam 
War. There are nearly 14,800 profiles (maximum number of observations in a year) in 1968 
(Figure 3.4). The minimum number of observations was in 1984 (near 3,00 profiles). There 
figures indicate that the number of observations within a given year exhibits a seasonal 
variation. 
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The main limitation of the MOODS data set is its irregular distribution in time and 
space. Certain periods and areas are over sampled while others lack enough observations to 
gain meaningful insight. Vertical resolution and data quality are also highly variable 
depending a great deal on instrument type and sampling expertise. 

We binned seventeen years’ temperature profiles (from 1968 to 1984) into monthly 
data for an individual year which produced 204 temperature data sets. For the salinity 
profiles, we binned into monthly climatology data bins which produced 12 mean monthly 
salinity data sets. 
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MOODS temperature profiles with 189,059 data 
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Figure 3.1. The Spatial distribution of MOODS temperature profiles 
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MOODS salinity profiles with 14,004 data 
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Figure 3.2 The Spatial distribution of MOODS salinity profiles 
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MOODS data temporal distribution Depth 0 M 
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Figure 3.3 The temporal distribution of MOODS 
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Figure 3.4 The temporal distribution of MOODS in 1968 
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IV. ESTABLISHMENT OF GRIDDED DATA 



For each temperature (or salinity) profile, we linearly interpolated the data into 
thirteen standard depths (Table 4. 1). After this process, we used the minimum curvature 
spline method [Smith and Wessel, 1990] from the Generic Mapping Tools (GMT) [Wessel 
and Smith, 1992] to create horizontal gridded data. This method interpolates the data (J>jj 
at observational points (xj, yj) onto the grid points with a spline function 4>(x,y) such that 
at observational points; 

,yj) = (^) 

The selection of the spline function is based on the minimum curvature principle, 
which requires continuous second derivatives and minimal total squared curvature, 

C = ff (V^(j>)' dx dy (1) 

Briggs [1974] showed that minimizing C is equivalent to solving the differential equation, 

V^(V2(J)) = yy-yj) (2) 

i 

where must be chosen such that Equation (1) should be satisfied and where 6 is the 
Kronecker delta function. To solve Equation (2) we need the following boundary conditions 
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( 3 ) 



dn^ 



= 0 



and 



A(v2({)) = 0 (4) 

on 



along the edges, where didn indicates a derivative normal to an edge, and 

= 0 (5) 

dx dy 



at the comers. Equations (3),(4),(5) are called free edge conditions; and with these 
conditions, equation (2) has a unique solution with continuous second derivatives called the 
natural bi-cubic spline. The nomenclature comes from an analogy with elastic-plate flexure. 
Small displacements (j) of a thin elastic plate of constant flexural rigidity D, subject to a 
vertical normal stress q and constant horizontal forces per unit vertical length of T^x, T^y, T^y, 
approximately satisfy 



i)V2(V2<j)) - [J 



a"(j) 



2T 



xy 



a^(j) 



=9 

’ dy'^ 



(6; 



[Love. 1927]. The minimum-curvature gridding equation (2) is a special case of equation 
(6) when horizontal forces are zero, and the boundary conditions represent zero bending 
moment on the edges (3), zero vertical shear stress on the edges (4), and zero twisting 
moment at the comers (5) [Timoshenko and Woinowsky-Krieger, 1968]. 
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Suppose that Txx= Tyy = T and T^y = 0, then Equation (6) becomes 



Z)V^(V^4>) - 7V^<t) = q 



( 7 ) 



When T = 0, Equation (7) is equivalent to Equation (2); but for an arbitrarily large T, the 
solution is dominated by the second term. Here, T has units of force per unit length and the 
T required to adjust the solution scales with D and q; Smith and Wessel [ 1 990] avoided this 
problem by writing 



where T, is a tension parameter and the I subscript indicates internal tension. When Tj = 0, 
equation (8) reduces to equation (2); and therefore, the minimum-curvature solution is one 
end-member case of equation (8). When 7) = 1, the first term in equation (8) vanishes; and 
the solution is harmonic between constraining points. One may prefer to think of this end 
member as representing the steady-state temperature field in a conducting plate with heat 
sources or sinks at the data points. For any 7} in 0< T;< 1, equation (8) gives a solution with 
continuous curvature, although it does not minimize equation (1) except when 7) = 0 [Smith 
and Wessel, 1990]. 

For the boundary condition, they implemented conditions (4) and (5) but replace 
condition (3) by 



(i-r,)v"(v"(|)) - r,v'(|> = , y-yj) 



( 8 ) 




+ = 0 



dn 



( 9 ) 
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where Tg is a tension parameter for the boundary which also varies between 0 and 1 . The 
free-edge condition corresponds to Tg= 0; Tg = 1 forces the solution to flatten at the edge 
(Figure 4.1a,b). This flattening is sometimes desirable, as when gridding potential anomalies 
which should decay toward a regional background field away fi-om the source region. 
Gridding with tension is an improvement because it adds a degree of fi'eedom and relaxes the 
minimum-curvature constraint. 
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a) 




Figure 4.1a Cross-sections through surfaces produced with splines 
in tension. The black squares are data constraints. The heavy line 
is the minimum-curvature end member, the thin line is the 
harmonic end member, and the dashed line is an intermediate 
case using some tension (from Smith and Wessel, 1990) 
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b) 
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Figure 4.1b Cross-sections through continental shelf and slope. The black 
squares represent the intersection between the measured bathymetry (dashed 
line) and 100m isobath contours. These intersections (contour coordinates) were 
then gridded using minimum curvature (heavy solid line) and some tension (thin 
solid line) The minimum-curvature method introduces an extraneous inflection 
point and exceeds the -100m level (The bathymetry in this region is bounded by 
the -200m and -100m level). The surface produced with tension gives a much 
better approximation since it suppresses local maxima and minima between data 
. constraints (from Smith and Wessel, 1990) 
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V. COMPOSITE ANALYSIS 



Our goal is to investigate the temporal and spatial variabilities of the SCS thermohaline 
structure. We use T,,iJ to represent temperature and salinity gridded data, where 

(Xj ,yj) are the horizontal grids, z* (= 1, 2,..., 13), is the vertical level, T, 
(=1968,1969,... 1984) is the time sequence in years, and (= 1, 2,..., 12) is the time 
sequence in month. We define the following temporal average; 

At I 



which is the long-term mean value for the month and 
= 1 - 

which is the ensemble (climatological armual) mean value. Here, at =17 years. 

A. ANNUAL MEAN 

The annual mean (1968-1984) temperature j field over the SCS shows the pattern 

of northeast-southwest oriented isotherms at the upper layer fi-om the surface to 75 m depth 
(Figure 5.1a). The armual surface mean has a rather weak horizontal temperature gradient, 
decreasing fi-om 28.5"C near the Borneo coast to 25°C near the southeast China coast. A 
stron temperature fi-ont is found near the Luzon Strait in the sub- surface layer fi-om 100 m 
to 400 m (Figure 5.1b). The SCS water temperature is quite uniform (e.g., near 15°C at 200 
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m) and much lower than the water east of the Luzon Strait, the West Pacific Ocean Water. 
This might imply that the annual mean Kuroshio intrusion is not significant in the subsurface 
layer. 

In the central SCS region at 15°N zonal cross-section (Figure 5.2a), the isotherms 
of T bend slightly downward from west to east and no frontal structure is found. However, 

in the northern SCS region at 20°N zonal cross-section, the isotherms of f show a strong 

frontal structure near 120° to 125°E (Figure 5.2b), which may suggest two distinguished 
water masses, the SCS water mass (both surface and under-surface^ and the West Pacific 
Ocean water mass. The front between the two water masses strengthens from surface to 400 
m depth. 

The aimual mean salinity S shows more complexity than temperature. In the upper 
layer (Figure 5.3a), a large amount of freshwater enters the SCS from the Zhujiang (Pearl) 
River in the northwest and from the Mekong River in the southwest. The Kuroshio brings 
the salty water through the Luzon Strait into the north SCS and forms a salty tongue (34 psu) 
stretching into the southeast China coast. As the depth increases, both the Kuroshio intrusion 
and the river run-off effects become less important. At 75 m depth, the Kuroshio intrusion 
effect is still evident, but not the river run-off effect. Below 200 m depth (Figure 5.3b), the 
salinity is more uniform (near 34.6 psu at 200 m) throughout the whole SCS. 

In the central SCS region at 15°N zonal cross-section (Figure 5.4a), the isohalines 
of 5 are almost parallel to the isobaths in the upper 100 m. Below 100 m depth, the salinity 
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is very uniform (34.4 psu). However, in the northern SCS region at 20°N zonal cross-section 
(Figure 5.4b), the 34.6 psu isohaline clearly shows the Kuroshio intrusion from the Luzon 
Strait into the northern SCS in the layer between 100 and 220 m reaching 1 15°E. 

B. SEASONAL VARIATION 

The long-term monthly mean values relative to the annual mean 

represent the composite features of the monthly mean variations of temperature and salinity. 

1. Temperature 
CL Sutface 

We choose the SST to represent the thermal feature of the surface layer. ' The 
seasonal SST variation obtained from the MOODS data is quite consistent with earlier 
investigations based on the National Center for Environmental Predictions (NCEP) data [Chu 
et al., 1997b]. 

During winter (January), T is negative almost everywhere throughout the 

whole SCS (Figure 5.5a). The typical winter (January) pattern contains 
northeast-to-southwest oriented isotherms in the northern SCS (north of 12° N), a warm 
anomaly (T > -1 .2°C) in the west of Borneo- Palawan islands (WBP), and a cool anomaly 
( r < -2.4°C) near the South Vietnam coast (SVC) in the southern SCS (south of 12°N). The 
spring-to-summer transition (March to May) pattern shows the northward expansion of the 
WBP warm anomaly and the formation of a central SCS (10°-15° N) warm anomaly with T 
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>0.6°C. The isoline 1°C encloses almost all the SCS in May. 



The summer (July) pattern is characterized by northeast-to-southwest oriented 
isotherms in the northern SCS (north of 15°N) (Figure 5.5b), a WBP cool anomaly ( f < 
0.4°C), and a SVC warm anomaly (> 1.2°C) in the southern SCS. The summer pattern is 
opposite to the winter pattern. The fall-to-winter transition (October) pattern shows the 
northward expansion of a WBP cool anomaly and the formation of a central SCS (10° - 
15°N) cool anomaly. 

b. Sub-Surface Level (75 m) 

We chose 75 m depth to represent uie sub-surface (Figure 5.6). The seasonal 
variation of this layer is characterized by the west Luzon coast (WLC) and SVC cool (warm) 
anomalies and southeast China coast (SCC) warm (cool) anomaly in winter (summer). The 

WLC cool anomaly appears in December, reaching its maximum strength in January with T 
< -3.4°C (Figure 5.6a), lasting two months and disappears in April. The WLC warm 

anomaly forms in June with T > 2°C, reaches T > 2.4°C in July (Figure 5.6b), and lasts for 

four months. 

Figure 5.7 shows the zonal cross-section at 19°N of thef field, representing ' 
the vertical structure of the WLC cool eddy in winter (Figure 5.7a), and the vertical structure 
of the WLC warm eddy in summer (Figure 5.7b). Both WLC winter cool and summer warm 
eddies reach 250 m depth. The winter cool-core ( T< -3 °C) and the sv -nmer warm-core ( T> 
2.5°C ) are located between 50 m to 100 m depth. 
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c. Intermediate Level (200 m) 

At 200 m depth, the seasonal variation of temperature is featured by a basin- 
wide cool pool (7’< -0.8°C), a strong SCC warm pool (f> 0.8°C), and a weak WBP warm 
pool ( r > 0.2°C) during the winter (Figure 5.8a) and a basin-wide warm pool ( 7 > 0. 4°C), 
and a weak SVC cool pool ( 7 < 0°C) during the summer (Figure 5.8b). The WLC winter 
cool eddy (7< -1.0°C ) and summer warm eddy (7>0.8°C ) are still quite evident at that 
depth. The weak WBP warm pool strengthens, expands toward the northeast in spring 
(Figure 5.8a), gradually occupies the whole SCS basin, and turns to the summer pattern in 
July (Figure 5.8b). In autumn, the weak SVC cool pool strengthens and expands toward 
northeast, gradually occupies the whole SCS basin and tums.to the winter pattern in January. 
Below 300 m depth, there is no significant seasonal variability. 

2. Salinity 

The seasonal salinity variation S is evident only in the layer above 200 m depth. 
Below this depth, the seasonal variability is very weak. 
a. Surface 

During winter (January), S is positive almost everywhere throughout the 

whole SCS except the SVC and the Gulf of Tonkin (Figure 5.9a). A salty tongue with.? >0.4 
psu stretches fi’om the Luzon Strait to the continental shelf along the southeast China coast, 
representing a strong Kuroshio intrusion. Another weak fi'esh pool with5 < -0.2 psu is found 
to the west of Luzon Island. During summer (July), a salty tongue with S around 0.4 psu 
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stretches from the Vietnam coast to 1 18°E between 10°-17°N, occupying the central SCS 
(Figure 5.9b). Everywhere else is occupied by negative S. A strong fresh tongue occupies 
the continental shelf along the southeast China coast and the Gulf of Tonkin. The negative 
values in the vicinity of the Luzon Strait indicate the cease of the Kuroshio intrusion into the 
SCS. On the other hand, the negative S in the southern SCS shows the influx of fresh water 
from the southern boundary. 

b. Sub-Surface Level (75 m) 

The seasonal 5 variation is featured by a wide salty tongue stretching 
southwestward from the Luzon Strait to the SCS basin in the wanter (Figure 5.10a) and a 
fresh tongue stretching northward from the southwest boundary near the Mekong River 
mouth to the SCS basin in the summer (Figure 5.10b). The former represents the intrusion 
of the salty Kuroshio water (with maximum <S > 0.6 psu) through the Luzon Strait, and the 
latter indicates the low salinity water (with minimum S < -0.4 psu) extending to the SCS. 

Figure 5.1 1 shows the zonal cross-section at 19°N of the S field, representing 
the vertical structure of the WLC high salinity anomaly (S> 0.4 psu) in wanter (Figure 5. 1 la) 
and low salinity anomaly {S< -0.2 psu) in summer (Figure 5.11b). Both WLC high /low 
salinity anomalies are located between 50 m to 100 m depth. 

c. Intermediate Level (200 m) 

At 200 m depth, the seasonal variation of 5 shows a similar pattern (Figure 
5.12) to that at 75 m depth. In the winter (Figure 5. 12a), positive S with a maximum value 
of 0.3 psu occupies the most area of the SCS, except the southwest comer near the Mekong 
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River, where S is negative (-0.4 psu). In the summer (Figure 5.12b), a narrow fresh tongue 
with a minimum value of -0.4 psu stretches northeastward from the southwest comer near the 
Mekong River to 20°N (Figure 5.12b). 
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Figure 5.1a The annual mean temperature field (f) 
from surface to 100 meters 
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Figure 5.1b The annual mean temperature field ( j) 
from 125 to 400 meters 
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Figure 5.2a The zonal cross-section view of annual mean temperature f 
from Latitude 14®N to Latitude 19®N 
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Figure 5.2b The zonal cross-section view of annual mean temperature j 
from Latitude 20®N to Latitude 25®N 
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Figure 5,3a The annual mean salinity field (s) 
from surface to 100 meters 
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Figure 5.3b The annual mean salinity field (s) 
from 125 to 400 meters 
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Figure 5.4a The zonal cross-section view of annual mean salinity ( 5 ') 
from Latitude 14°N to Latitude 19°N 
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Figure 5,4b The zonal cross-section view of annual mean salinity (s) 
from Latitude 20°N to Latitude 25®N 
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Figure 5.5a Monthly mean temperature anomalies ( at surface 

from January to June 
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Figure S.5b Monthly mean temperature anomalies ( at surface 

from July to December 
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Figure 5.6a Monthly mean temperature anomalies ( T) at 75 meters 

from January to June 
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Figure 5.6b Monthly mean temperature anomalies ( T) at 75 meters 

from July to December 
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Figure 5.7a Zonal cross-section view of monthly mean temperature 
anomalies ( 7’) at 19®N from January to June 
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Figure 5.7b Zonal cross-section view of monthly mean temperature 
anomalies ( T) at 19°N from July to December 
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Figure 5.8a Monthly mean temperature anomalies ( f) at 200 meters 

from January to June 
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Figure 5.8b Monthly mean temperature anomalies ( T) at 200 meters 

from July to December 
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Figure 5.9a Monthly mean salinity anomalies (S) at surface 
from January to June 
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Figure 5.9b Monthly mean salinity anomalies (S) at surface 
from July to December 
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Figure 5.10a Monthly mean salinity anomalies (S) at 75 meters 
from January to June 
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Figure 5.10b Monthly mean salinity anomalies (S) at 75 meters 
from July to December 
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Figure 5.11a Zonal cross-section view of monthly mean salinity 
anomalies (S) at 19°N from January to June 
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Figure 5.11b Zonal cross-section view of monthly mean salinity 
anomalies (S) at 19°N from July to December 
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Figure 5.12a Monthly mean salinity anomalies (S) at 200 meters 

from January to June 
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Figure 5.12b Monthly mean salinity anomalies (S) at 200 meters 
from July to December 
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VI. EMPIRICAL ORTHOGONAL FUNCTION (EOF) ANALYSIS 



The empirical orthogonal function (EOF) analysis is a useful tool to describe the main 
mode of temporal and spatial variabilities. It is usually applied to the anomaly data relative 
to the climatological mean. 

A. THERMAL VARIABILITY 

The temperature anomaly field relative to the climatological annual mean field 

- f{x^,y.,z^) (13) 

is re-arranged into a NxP matrix , /~ ) , n = 1, 2, ..., N ; and p =1, 2 ,..., P. Here P 

= 204 is the total number of time points used for computing the covariance matrix, that is, 1 7 
years of monthly data sets; 7/ = 5733, corresponds to the number of grids (i = 1, 2, ..., 21;y 
= 1, 2,..., 21; A:= 1,2,..., 13). The non-seasonal thermal variability is investigated using the 
Empirical Orthogonal Function (EOF) analysis. This method separates the data sets into 
eigenmodes. Generally speaking, each mode has an associated variance, a dimensional spatial 
pattern, and a non-dimensional time series. For example, Lorenz (1965) used the eigenvectors 
of the covariance matrix associated with a sample data field to calculate the time series of 
coeflBcients or predictors for the statistical weather prediction. Chu et al. (1997a, b) used this 
method to investigate the non-seasonal variabilities of SST and surface wind stress. From 
our data , a 5733x5733 spatial covariance matrix is calculated by 
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^11 ^12 ••• 



R = 



^21 ^22 



2N 



, R = 

~ nrn 






^TnrJ^) t(rj^), A'=5733, P=204 (14) 

^ P 



^N1 ^N2 ••• ^NN 

where n and m (1,2, ...,5733) denote the grid locations. The diagonal elements of the 
covariance matrix =(n=l,2,...,N) are the variance at location r„ The ofif-diagonal elements 
are the covariance with spatial lag equal to the difference between the row and column 
indices. This symmetric matrix has its 5733 real eigenvalues , and eigenvectors (rj) 
such that 






(15) 



The eigenvectors 3>j, Oj, ..., are called Empirical Orthogonal Functions with the unit of 
°C. Each is a 5733-point (21x21x13) distribution of non-seasonal variation temperature 
anomaly pattern. The eigenvalues, X„(a=l, 2, ..., N), are all positive and the summation of 
them, EXjj, equal the total variance. Therefore, is considered as portion of total variance 
"explained" by the EOF . It is convenient to label the eigenfimctions, so that the 
eigenvalues are in descending order, i.e., Xj > > Xj > .... 

The data matrix f(r^ , t^) is thus written approximately by 
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where PCJZi, tj is the principal component representing the temporal variation of the 
associated spatial pattern described by the EOF mode, 0„(x,, yj,zJ. 

B. PRINCIPAL EOF MODES 

The first six leading EOF modes are able to account for 64. 1% of the total variance 
(Table 6.1). The first two leading modes, EOFl and EOF2, explain 44.4% of the total 
variance. Each EOF mode is three-dimensional and normalized so that its total spatial 
variance is equal to unity. We may easily see the vertically coherent structure of each mode. 

1. The EOFl Mode 

The EOFl mode (<?i) shows frontal and multi-eddy structures. The frontal structure 
is evident in the northern SCS. Only positive values occur in the upper layer (above 75 m 
depth) and both positive and negative values appear at the other depths (Figure 6.1a,b). 
However, we cannot use these positive/negative values of 0i to identify the warm or cool 
anomaly of the individual eddy. This is because the temperature anomaly, due to the EOFl 
pattern, is the product of PCj and 0i. 

a. Western Boundary Current and Frontal Structures 
From the surface to 30 m depth, a bi-frontal structure (northeast-to-southwest 
oriented isolines) is the major thermal feature in the northern SCS (north of 12° N): a 
strong coastal front along the China coast, and a relatively weak and wide front across the 
SCS basin from the Vietnam coast to the Luzon island between 12°-17°N. The two fronts 
change their strengths synchronously, since the values of 0^ are positive for the two fronts. 
The bi-frontal structure weakens with the depth and disappears at 50 m depth. The front 
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reoccurs from 75 m depth and continues to be seen down to 300 m depth. Below 75 m 
depth in the northern SCS, shows a northeast-to-southwest oriented single front structure 
aligning to the China coast. The values of across the front change signs from negative 
near the China coast to positive in the middle of the northern SCS, indicating the out-of- 
phase variability of the front near the China coast and the deep basin. 

Furthermore, we may notice that the values of <?j change signs from the top 
layer (0-30 m) to the lower layer (100-300 m) near the western boundary. This implies a 
strong baroclinicity of the western boundary currents with the opposite directions at the top 
layer and the lower layer. 

b. Southern SCS Multi-Eddy Structure 

In the southern SCS (south of 12°N) from the surface to 30 m depth, the 
EOFl shows a dipole structure (dual eddies) with a high center (0.3°C) near the mouth of 
the Mekong River, and a low center (< 0. 14°C) in the west of the Bomeo-Palawan islands 
(WBP). This dipole structure is still quite evident from 30 m to 75 m depths. Below 75 m 
depth, the western eddy, near the Mekong River mouth, is much stronger that the eastern 
eddy. 

At 30 m depth, an eddy (>0.4°C) appears in the northern SCS. This eddy 
slants towards the east (the Luzon Strait) with depth and becomes very strong in the layers 
between 75 m and 150 m. Below 150 m depth, this eddy reduces its strength. 
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c. 



Kuroshio Intrusion 



In the northern SCS, the isolines of 0^ intersect the cross-section of the 
Luzon Strait with the patterns showing Kuroshio intrusion/SCS outflow through the Luzon 
Strait. This pattern indicates a weak Kuroshio intrusion/SCS outflow in the upper layer 
above 50 m depth, becomdng stronger below 50 m depth and reaching maximum strength 
in the layer between 100m and 150 m depths. 

2. The Second EOF Mode 

Similar to the first EOF mode, the second EOF mode {0^ has only positive values 
in the upper layer from the surface to 50 m depth and has both positive and negative values 
appearing on the other depths. The negative 0j area is limited in the southeast comer of the 
SCS with small absolute values (Figure 6.2a,b). This means that the overall SCS basin 
thermal variability is either warming or cooling. We may call this mode the 
warming/cooling mode. In addition, the isolines of 0j at all the depths near the Luzon Strait 
are almost parallel to the cross-section of the strait. This indicates that the second EOF 
(warming/cooling) mode is not caused by the Kuroshio intrusion/SCS outflow. It might be 
related to the atmospheric forcing. Furthermore, the second EOF mode also shows coastal 
front and multi-eddy stmctures. 

a. Coastal Front Structure 

In the north SCS, a coastal front structure can be identified on the continental 
shelf near the China coast (Figure 6.2a). As the depth increases, this front stmcture extends 
to the southern SCS and occupies the western part of the SCS from 100m to 300 m depths 
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(Figure 6.2b). This strong gradient of 02 caused by the western boundary current. 

b. Northern SCS Multi-Eddy Structure 

The EOF2 clearly shows a multi-eddy structure in the northern SCS sub- 
surface layer between 50 m and 125 m. At 50 m depth, a strong north-south dipole structure 
(dual eddies) with a high center (0.2°C) in the west of the Luzon Strait (18°-22°N, 1 15°- 
1 19° E), and a low center (< 0.04°C) in the central SCS (12°-15°N, 1 10°-1 16° E) . This 
dipole structure strengthens at 75 m depth and weakens as depth increases from 75 m. 
Below 125 m depth, the northern SCS multi-eddy structure disappears. 
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e 6.1 Variances of the First Six Leading EOFs 



C. TEMPORAL VARIABILITIES 

Each principal component can be treated as the projection of the temperature 
anomaly through a ‘filter’ of corresponding EOF mode (Chu et al., 1997). 
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1 . 



First Principal Component 



The first principal component, , is the time series (1968-1984) of 

the EOFl amplitudes and shows basically a seasonal variation. During the winter (summer) 
monsoon season, has negative (positive) values (Figure 6.3), which indicate 

that the temperature anomaly, due to this mode, has the opposite (same) sign with respect 
to EOFl . Since the first EOF mode is always positive throughout the upper layer (0 - 
50 m) of the whole SCS (Figure 6.1), positive (negative) values of PC,(/~) during 

summer (winter) monsoon seasion correspond to warm (cool) anomalies. This shows the 
regular seasonal variation.’ 

An interesting feature of the first EOF mode 0^ is its n^ative value near the western 
boundary in the sub-surface layer (75 - 300 m). Since F*Cj (^) has negative (positive) 

values during winter (summer) monsoon season, the thermal variability of the western 
boundary is out-of-phase with respect to the seasonal variability, that is, warm (cool) 
anomaly in the winter (summer) monsoon season. The warm anomaly in the winter could 
be related to the Kuroshio intrusion. 

Accordingly, the seasonal variation represented by PC^it^) is not a pure 

sinusoidal. The summer had a larger variability in the 1980's than in the early 1970's. 
However, the winter had a larger variability during 1968-1972 than the rest of the periods. 
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2. Second Principal Component 

The second principal component, PC^ ( ) , is the time series (1968-1 984) of the 

second EOF (warming/cooling) mode amplitudes and shows interannual variabilities with 
periods of two to five years (Figure 6.3). The positive (negative) values v~f 
PCj indicate the SCS warming (cooling). The maximum value of 02 is 0.3°C near 

the China coast at the surface and a high center with isoline of 0.2°C in the west of Luzon 
Strait at depths between 50 m aid 125 m (Figure 6.2). 

The values of PC^itp) between 10 and -10. Its maximum value (10) 

appeared in August 1969 and September 1976, which indicates the strong SCS warm 
anomaly with the maximum warming of 3°C [0.3°Cx 10] at the surface near the China coast, 
and the occurrence of a sub-surface (50-125 m) warm eddy with the anomaly of 2°C 
[0.2°Cx 10] in the west of the Luzon Strait. Its minimum value (-10) appeared in October 
1974, which indicates the strong SCS cool anomaly with the maximum cooling of -3°C 
[0.3°Cx(-10)] at the surface near the China coast, and the occurrence of a sub-surface (50- 
125 m) cool eddy with the anomaly of -2°C [0.2°Cx(-10)] in the west of the Luzon Strait. 

We foimd the following criterions in evident warm/cool anomaly occurrence with the 
maximum temperature anomaly exceeding 1.5°C. When PC 2 (tp) ^ the SCS has an 

evident warming (maximum 1.5°C) and a warm eddy occurs at the west of the Luzon Strait 
in the sub-surface with the warm anomaly of 1 °C. On the other hand, />Q ( ^ ) < 5 , the 
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SCS has an evident cooling (maximum above 1.5°C) and a cool eddy occurs at the west of 
Luzon Strait in the sub-surface with the cool anomaly of -1°C. 

The strong anomaly periods (maximum above 3°C) and evident periods (maximum 
above 1.5°C) are listed in Table 6.2. We may find that the SCS warm anomaly generally 
appeared during the period of 1968-1972. After 1972, the cool anomaly occured more often, 
except in September of 1976. 



Anomaly 

Period 


Maximum Anomaly 
Above 3®C 


Maximum Anomaly 
Above 1.5®C 


SCS Warming 


Aug. 1969 
Sep. 1976 


Sep. 1968, 

Jan. Feb. May Jul. Aug. 1969 
Nov. 1969 
Jun. 1970 
Oct. Dec. 1971 
May Jul. Sep. Nov. 1972 
Sep. 1976 
Jul. 1977 


SCS Cooling 


Oct. 1974, 


Oct. 1974 
Jan. 1976 
Feb.- Mar. 1977 
Feb. 1980 
Apr. 1983 



Table 6.2 SCS warm/cool anomaly periods 
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Figure 6.1a EOF 1 Temperature anomalies (Unit 0.01°C) 
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Figure 6.1b EOF 1 Temperature anomalies (Unit 0.01°C) 
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Figure 6.2a EOF 2 Temperature anomalies (Unit 0.01°C) 
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Figure 6.2b EOF 2 Temperature anomalies (Unit 0.01°C) 
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Figure 6.3 The first six principal components (PC) for 1968-1984 
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vn. THE SOUTH CHINA SEA CmCULATION 



A. GEOSTROPmC VELOCITY INVERTED FROM DENSITY FIELD 

The thermal wind relation is often employed to calculate geostrophic velocity fi'om 
hydrographic data 




where (u, v) and (u^, v^) are the geostrophic velocity at any depth z and at a reference depth 
Zo, p is the in-situ water density, Po is the characteristic value of the density, f=2QsinA 
is the Coriolis parameter, where Q is the earth's rotation rate and X is the latitude. The 
hydrographic data determine the geostrophic vertical shear. The reference velocity (u^, vj 
can be determined by several inverse techniques, such as the P-spiral method [Stommel and 
Schott, 1977; Schott and Stommel, 1978], the box inverse method [Wunsch, 1978], and the 
Bernoulli method [Killworth, 1986]. Wunsch [1994] proposed a generalized p-spiral and 
geostrophic balance inverse method based on a nonlinear optimization on Needler's P 
equation [Needier, 1967]. Davis [1978] pointed out that the P-spiral method and the box 
inverse method, regardless of appearances, are based on the same dynamical assumptions, 
and the existence of the P-spiral thus becomes a required condition for these inverse 
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methods. 



Non-geo strophic inverse methods were also developed. Bogden et al. [1993] 



estimated the time-averaged velocity field in the North Atlantic fi'om observations of density, 
wind stress and bottom topography based on the assumption of minimum mixing of potential 
density at mid-depth. Tziperman et al. [1992] applied the adjoint method to investigate the 
steady-state oceanic general circulation using a realistic numerical oceanic general circulation 
model and hydrographic data. 

Recently, Chu [1995] and Chu et al. [1998a] proposed a simple P-vector inverse 
method to obtain the absolute velocity from hydrographic data. A velocity vector is 
represented into the product of a unit vector (P-vector) and a scalar. 



The inversion of V is fulfilled by two steps: (1) determination of the unit vector P, and 
(2) determination of the scalar r(x,y,z). Since the geostrophic flow is along both the 
isopycnal and the iso-potential vorticity surfaces, the velocity V is parallel to Vq x Vp. Here, 
p is the potential density, and 



is the potential vorticity. If the p surface is not parallel to the q surface, the unit vector P cc . 
be defined by [Chu, 1995] 



V = r(x, y, z) P. 



(19) 



q = f dp / dz 



( 20 ) 



p ^ Vp X Vg 

|Vp x V^l 



( 21 ) 
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Applying the thermal wind relation at two different depths and z^, as shown in 



Figure 7. 1, a set of algebraic equations for determining the parameter r is obtained 



where p is the in-situ density. 

Chu et al. [1998b] shows the advantages of using the P-vector method, namely, the 
easy identification of the two necessary conditions for the P-spiral inverse method: (1) 
intersection between the p- and q-suifaces (existence of the vector P), and (2) the P turning 
(turning of the vector P). 

For water columns satisfying the two necessary conditions, equation (22) may be 
used to compute for the level z^. There are (N-1) sets (m = 1, 2, k-1, k+1, ... N) of 
equations (22) for calculating r^l Here, N is the total number of vertical levels of the water 
column. All of the (N-1) sets of equations are compatible under the thermal wind constraint 
and should provide the same solution. However, due to errors in measurements 
(instrumentation errors) and computations (truncation errors), the parameters r^^ may vary 
with m. A least square error algorithm [Chu et al., 1998b] is used to minimize the error. 
For further details and validation of the algorithm presented, see Chu et al. [1998a]. 




( 22 ) 



which are two linear algebraic equations for r^^ and r^”^ [r ® = r(x,y,Zj )]. Here 




(23) 
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B. SEASONAL VARIATION OF SCS CIRCULATION 



1. General Description 

The most obvious features of the inverted summer and winter SCS circulation are 
flow separation, basin gyres, and mesoscale eddies. Hinted in Wyrtki's (1961) depiction 
but more explicit in our computation, are the upper level (10 m depth) west-to-east 
cross-basin currents, which are almost parallel to Ib'N latitude fi-om February to July. The 
inverted monthly SCS circulation agrees quite well with earlier observational study (Wyrtki, 
1961). As the commonly used inverse methods, theP-vector inverted western boundary 
currents are somewhat small. This could be caused by the neglect of relative vorticity in the 
potential vorticity calculation [see (21)]. 

2. Surface Circulation 

During the winter monsoon period (November to March), the winter Asian high 
pressure system brings strong winds from the northeast (Figure 2.4) and the SCS surface 
circulation pattern is cyclonic (Figure 2.5). At the surface inflow from the Luzon Strait (the 
Kuroshio intrusion) turns southwest along the Asian continental shelf, then south along the 
coast of \^etnam (Figure 7.2). In November, this current splits into two currents at 12°N: 
a southward alongshore current and an eastward off-shore current. The southward branch 
turns east at 6°N, then turns northeast and joins the eastern branch at 1 15®E. An evident 
multi-eddy structure is obtained with two cyclonic eddies in the north SCS (north of 12°N) 
and cyclonic and anticyclonic eddies in the southern SCS (south of 12°N). The flow pattern 
does not change vary much from November to December, except the cyclonic eddy near the 
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Luzon Strait strengths. In January, the west cyclonic eddy in the northern SCS expands 
toward the south. In February, the west cyclonic eddy in the northern SCS disappears and 
instead a cross-basin flow appears along 16°N latitude and joins the cyclonic eddy near the 
Luzon Strait. The southern SCS is occupied by an anticyclonic eddy. This anticyclonic eddy 
expands to north and forms the only evident eddy in the central SCS in March (Figure 7.2a). 
This eddy is associated with the central SCS warm pool in spring (Appendix A) and expands 
toward the south in April. May is the month of the summer monsoon transition. The 
anticyclonic eddy reduces its strength and becomes unidentifiable. 

During the summer monsoon period (Mid-May to August), winds blow fi'om the 
southwest (Figure 2.2) and the SCS surface circulation generally follows suit with 
anticyclonicity in the southern basin (Figure 2.3). Inflow is through the southern boundary 

f 

and outflow is through the northern Taiwan Strait and eastern Luzon Strait. The inverted 
July means the general circulation pattern has the following features. Velocities reach 20 
cm/s near the Vietnam Coast and splits into two currents at 12°N: the coastal current and 
ofi-shore current. The off-shore current fiarther bifurcates and partially leaves the coast; the 
bifurcation point is at 1 10°E. The northern branch moves northeast to 1 13°E and then turns 
zonally between 15°-18°N. The southern branch moves zonally until reaching the Palawan 
Island, then turns north and joins with the north branch at 16°N. A cyclonic eddy appears in 
July near the south Vietnam coast and strengthens in August. 
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3. 



Sub-Surface (75 m) 



The sub-surface (75 m) shows a strong seasonal variation (Figure 7.3). Warm and 
salty Kuroshio water enters the SCS through the Luzon Strait in October-March, the 
transition to and during the winter monsoon season. The winter circulation pattern is the 
basin-wide cyclonic gyre (November - January). A cyclonic mesoscale eddy splits from the 
gyre near the Luzon Strait in January and keeps its strength in February. On the other hand, 
the basin- wide cyclonic gyre weakens and shrinks toward the southwest in February and 
disappears in March. 

During the summer monsoon period (June-September), the Kuroshio intrusion 
through the Luzon Strait ceases. A weak anticyclonic eddy occurs in June near the South 
^^etnam coast and strengthens in July. A cross-basin current establishes in August from the 
north branch of that eddy to the Luzon Strait and forms a weak cyclonic gyre in the 
northwest and an anticyclonic gyre in the southeast. During the winter monsoon transition 
period (October), the southeast ^ticyclonic gyre weakens and the northeast cyclonic gyre 
gains strength and eventually the cyclonic gyre occupies the whole SCS basin, the beginnmg 
of the winter circulation pattern. 

4. Intermediate Level (400 m) 

As the water depth increases, the inverted velocity generally reduces with depth to 
200 m. The velocity at 200 m depth is quite small (see Appendix D). However, the velocity 
increases with depth below 200 m. Figure 7.4 shows the monthly variation of the SCS 
circulation at 400 m. At that depth, the inverted flow flows out of the SCS through r e 
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Luzon Strait. The winter pattern is characterized by three anticyclonic eddies (December- 
February), with two eddies in the northern SCS, and one in the southern SCS. The 
anticyclonic eddy near the Luzon Strait is evident from December to February. An east-to- 
west cross-basin current establishes along 16°N in February from the southern branch of 
the Luzon eddy. This cross-basin current sustains from February to April and separates the 
SCS basin into two parts. The northern part is occupied by an anticyclonic eddy and the 
southern part by a cyclonic eddy. The summer circulation pattern is featured by double 
anticyclonic gyres. The southern gyre is quite evident in June and July. The tangential 
velocity is around 4 cm/s. 

5. Kuroshio Intrusion Through the Luzon Strait 

The SCS joins the Pacific Ocean through the Luzon Strait. The seasonal variations 
of the intrusion of the Kuroshio Water into the SCS through the Luzon Strait have been 
investigated in earlier studies [Fan and Yu, 1981; Shaw, .1989; Shaw, 1991]. As pointed out 
by Shaw [1991], the intrusion process is important not only to the transport of water masses 
into the SCS, but also to the shelf circulation off the southeast coast of China. Shaw [1989, 
1991] used the discriminant analysis method to classify the water mass T, S characteristics 
at 150, 200, and 250 meters, and found that water characteristics of the Philippine Sea 
(Kuroshio) was identified along the continental margin south of China from October to 
January. The presence of this water indicated an intrusion current from the Philippine Sea 
info the SCS. 

To verify Shaw’s result, we computed the east-west geostrophic velocity along 
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120.5°E longitude using the P-vector method. The temperature and salinity used for the 
computation are the gridded data from the surface to 400 m depth along 120° and 121°E 
(Figure 7.5). The monthly zonal velocity in 120.5°E cross-section is given in Figure 7.6. 
The negative values indicate westward flow through the Luzon Strait (Kuroshio intrusion). 
The positive values show the outflow from the SCS to the Philippine. Sea. Flow through the 
Luzon Strait has a strong seasonal variation, as well as a vertical shear. Evident Kuroshio 
intruaon (negative values) occurs during the winter monsoon season (November - March) 
in the upper 200 m. Below 200 m, the velocity is quite small, and most values are positive. 
We also computed the upper layer (0-400 m) volume transport through the Luzon Strait 
(Figure 7.7). The evident Kuroshio intrusion appears from November to March. 
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(a)Thermal wind relation at two different depths z, and z 
(b) p spiral turning angle (a^) between the two levels z, and z 
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Figure 7.2a Surface Circulation from Jar jary to June 
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Figure 7.2b Surface Circulation from July to December 
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Figure 7.3a Sub-Surface (75 m) Circulation from January to June 
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Figure 7.3b Sub-Surface (75 m) Circulation from July to December 
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Figure 7.4a Intermediate Level (400m) Circulation 
from January to June 
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Figure 7.4b Intermediate Level (400 m) Circulation 
from July to December 
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Figure 7.5 Temperature and Salinity data from surface to 400 meter used 
for the computation of east-west geostrophic velocity along 120.5®E 
(dash lines: T and S; solid line: geostrophic velocity) 
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Figure 7.6a The monthly zonal velocity in 120.5E from January to June 
( negative value indicate westward flow; Kuroshio intrusion) 
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Figure 7.6b The monthly zonal velocity in 120.5E from July to December 
( negative value indicate westward flow; Kuroshio intrusion) 
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Mass Transport (Sv) at the Luzon Strait from MOODS 




Figure 1.1 The Mass Transport (Sv) at the Luzon strait 
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vm. CONCLUSIONS 



The goals of this study were to investigate the thermohaline structure and the 
circulation patterns of the South China Sea using the historical MOODS hydrographic data. 
We applied the minimum curvature with spine method to produce thirteen level gridded data 
from the MOODS. Although the MOODS contain more than 1 80,000 temperature profiles, 
the deeper layer (below 400 m depth) is poorly sampled. Futhermore, the MOODS contain 
only 14,000 salinity profiles, so we were only able to establish climatological annual and 
monthly mean salinity fields at. thirteen levels above 400 m depth. After the establishment 
of the gridded temperature and salinity data, we used the composite and EOF analysis to get 
the seasonal and interannual variabilities of the SCS thermohaline structure, and the P-vector 
inverse method to obtain the three-dimensional absolute geostrophic velocity fields. The 
major results from this research are summarized as follows. 

A. CLIMATOLOGICAL MEAN THERMOHALINE FEATURES 

The annual mean (1968-1984) temperature f field over the SCS shows the pattern 
of northeast-southwest oriented isotherms at the upper layer from the surface to 20 m depth 
with a weak gradient from 28.5°C near the Borneo coast to 25°C near the southeast China 
coast. A strong temperature front is found near the Luzon Strait in the sub-surface layer 
from 100 m to 400 m. The SCS water temperature is quite uniform (e.g., near 13°C at 250 
m) and much lower than the water east of the Luzon Strait, the West Pacific Ocean Water. 
The aimual mean salinity ~S shows that in the upper layer, a large amounts of freshwater 
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enter the SCS from the Zhujiang (Pearl) River in the northwest and from the Mekong River 
in the southwest. The Kuroshio brings the high salinity water through the Luzon Strait into 
the north SCS and forms a salty tongue (34 psu) stretching into the southeast China coast. 
B. SEASONAL THERMOHALINE VARIABILITY 

Seasonal thermohaline variability is associated with fronts and multi-eddies. The bi- 
frontal structure is the major thermal feature in the northern SCS (north of 12°N): a strong 
coastal front along the China coast, and a relatively weak and wide front across the SCS 
basin from the Vietnam coast to the Luzon island between 1 2°- 1 7°N. The two fronts change 
thdr strengths synchronously; they weaken with the depth and disappear at 50 m depth. The 
front reoccurs from 75 m depth and is evident as deep as 300 m depth. At 30 m depth, an 
eddy (>0.4°C) appears in the northern SCS and slants toward the east (the Luzon Strait) 
with depth and becomes very strong in the layer between 75 m and 150 m. Below 150 m 
depth, this eddy reduces its strength. In the southern SCS (south of 12°N) from the surface 
to 30 m depth, a dipole structure (dual eddies) was found with a warm (cool) eddy near the 
Mekong River mouth, and a cool (warm) eddy in the west of the Borneo- Palawan islands 
during the summer (winter) monsoon season. Bdlow 75 m depth, the western eddy near the 
Mekong River mouth is much stronger that the eastern eddy. 

The seasonal salinity variability was also identified. During winter (January), a salty 

tongue 1 .2 psu anomaly stretches from the Luzon Strait to the continental shelf along the 
southeast China coast, representing a strong Kuroshio intrusion. A weak fresh anomaly (- 
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0.2 psu) is found in the west of the Luzon island. During the summer (July), a salty tongue 
with 0.4 psu stretches from the Vietnam coast to 1 18°E between 10°-17°N, occupying the 
central SCS while the rest is occupied by negative anomaly. A strong fresh tongue occupies 
the continental shelf along the southeast China coast and the Gulf of Tonkin. 

C KUROSmO INTRUSION 

Both water mass analysis and the inverted geostrophic velocity at the Luzon Strait 
confirm the Kuroshio intrusion during the winter monsoon season. Such an intrusion has a 
vertical shear. Evident Kuroshio intrusion (negative values) occurs during the winter 
monsoon season (November - March) in the upper 200 m while below 200 m the velocity 
is quite small and outflows from the SCS into the Philippine Sea. The computed upper layer 
(0-400 m) volume transport through the Luzon Strait has a strong seasonal variation; the 
maximum Kuroshio intrusion into the SCS with a value of -2.2 Sv occurring in March, and 
the maximum SCS outflow with a value of 1.5 Sv in September. 

D. INTERANNUAL SCS WARMING/COOLING PHASES 

The interannual SCS warming/cooling phases were identified by the second EOF 
mode with periods of two to five years. The strongest SCS warming occurred in August 
1969 and September 1976, with a warming of 3°C at the surface near the China coast and 
an occurrence of a sub-surface (50-125 m) warm eddy (2°C) in the west of Luzon Strait. 
The strongest SCS cooling appeared in October 1974 with a cooling of 3°C at the surface 
near the China coast, and an occurrence of a sub-surface (50-125 m) cool eddy (-2°C) in the 
west of the Luzon Strait. 
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The warming events appeared between the late 1960's and first half of the 1970's, and the 
cooling events occurred in the second half of the 1970's and early 1980's. 

E. MULTI-EDDY STRUCTURE 

The T, S and the inverted velocity fields show the multi-eddy structure: two cyclonic 
eddies in the north SCS (north of 12°N) and cyclonic and anticyclonic eddies in the southern 
SCS (south of 12°N)- The anticyclonic warm-core eddy expands to the north and forms 
the only evident eddy in the central SCS in March, and expands toward the south in April. 
May is the month of the summer monsoon transition. The anticyclonic warm-core eddy 
reduces its strength and becomes unidentifiable. During the summer monsoon period (mid- 
May to August), winds blow fi-om the southwest and the SCS surface circulation generally 
follows suit with anticyclonicity in the southern basin. Inflow is through the southern 
boundary and outflow is through the northern Taiwan Strait and eastern Luzon Strait. The 
inverted July mean general circulation pattern has the following features. Velocities reach 
20 cm/s near the Vietnam Coast and splits into two currents at 12T^: the coastal current and 
oflf-shore current. 
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